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ABSTRACT 

We propose a scenario to account for the surprising orbital properties of the 
planet orbiting the metal poor red horizontal branch star HIP 13044. The or- 
bital period of 16.2 days implies that the planet went through a common envelope 
phase inside the red giant branch (RGB) stellar progenitor of HIP 13044. The 
present properties of the star imply that the star maintained a substantial enve- 
lope mass of O.3M , raising the question of how the planet survived the common 
envelope before the envelope itself was lost? If such a planet enters the enve- 
lope of an RGB star, it is expected to spiral-in to the very inner region within 
< 100 yr, and be evaporated or destructed by the core. We speculate that the 
planet was engulfed by the star as a result of the core helium flash that caused 
this metal poor star to swell by a factor of ~ 3 — 4. The evolution following 
the core helium flash is very rapid, and some of the envelope is lost due to the 
interaction with the planet, and the rest of the envelope shrinks within about 
a hundred years. This is about equal to the spiraling-in time, and the planet 
survived. 



INTODUCTION 



In a recent paper ISetiawan et all ( 120101 ) announced the detection of a planet orbiting 
the metal-poor red horizontal branch star (HB) HIP 13044 (CD-36 1052) with an orbital 
period of P = 16.2 ± 0.3 days. The star resembles red HB stars in globular clusters, hav- 
i ng an effective temperature, a mass, a r adius, and a metalli city of T e ff = 6025 ± 63 K 



( ICarney et al. 



2008M IRoederer et aUBOlOh . M, 



6.7 ± (ICarney et all 2008b) , and [Fe /Hl 



0.8 ± (ISetiawan et aUBOlOh. R* 
-- -2.1 (feeers et all Il990l : Ichiba fc Beers 



2000l ; ICarney et al.l l2008bt IRoederer et al.l 120101 ). respectively. From the perspective of 
known exoplanets around main sequence stars, planets are rare around metal poor stars 
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[e.g., iSozzetti et al.l 120091 ). The detection of a planet around such a metal poor low mass 
star might be taken as a surprise, although for more than a decade theoretical studies have 



been proposing the pr e sence of planets in globular clusters (poker 



Soker fc Harpazl |2000L ISoker fc Harpazl 120071 : ISoker fc Hadar 



1998 



2001 



Siess fc Livio 



1999 



Soker fc Hershenhorn 



2003). 



What we find surprising are the orbital semimajor axis of a = 0.116 ± 0.01 AU and 
its eccentricity of e = 0.25 ± 0.05. These parameters raise the following questions, (i) How 
come a companion of a mass of only M p ~ 7.5Mj up survived a common envelope phase with 
a RGB star that did not lose its entire envelope? The properties of the star imply that the 
present envelope mass is M env ~ O.3M . Naively, one would expect that an envelope that is 
~ 40 times as massive as the planet would have caused the planet to continue spiraling-in 
inside the progenitor RGB envelope within a very short time (^3]). (ii) How come the orbit 
of a low mass companion that emerges from a common envelope is eccentric? Gravitational 
drag and tidal interaction with the envelope are expected to circularize the orbit. 

In section [2] we show that the angular momentum of the system is about equal to the 
initial (pre- RGB phase) orbital angular momentum of the planet. Therefore, there is no 
need to postulate the presence of a third body in the system. In section [3] we propose that 
the interaction of the planet with the RGB envelope took place over a relatively short time 
of about a hundred years. The interaction we speculate about in section [4] was triggered by 
a brief but substantial expansion of the star as a result of the core helium flash. In section 
H]we also summarize and conclude that some globular clusters should be a prime target for 
the search of planets around metal poor stars. 



ANGULAR MOMENTUM CONSIDERATIONS 



The possibility that planets spin-up RGB stars goes back to lPeterson et al.l (119831 ). who 
try to account for the fast rotation of s o me HB stars. New e r claims for planet-induce d 
RGB stellar rotation include Soker jl998h . lNelemans fc Tauri] Jl99gh . Isiess fc Livid ( jl99gh : 
Reddy et al.l J2OO2]), benissenkov & HerwisJ J2004h. iMassarottil J2008I) . ICarnev et all (£51, 
2008af ) (who i nclude HIP 13044). and ICarlberg et al.l (|2009l. boioh. A systematic study was 
conducted by lSoker fc Harpazl (120001 ) . whose calculations, assumptions, and approximations 
we adopt. 



The parameters we use here are as given and derived by ISetiawan et al.l ( 2010T). The 
compa nion mass is M p = 1.25 ±0.05Mj up / sinz, where Mj up is Jupiter mass. ISetiawan et al. 
( I2010l ) adopt the mean activity period to be due to rotation with P rot = 5.53±0.73 days, and 
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from that deduced the inclination angle of the orbital plane to be i = 9.7° ± 1.3. We note 
that, within the uncertainties, the ratio of orbital to rotation period equals 3. The possibility 
of a tidal resonance should be examined. We will scale qua ntities with % = 10°, a nd hence 
with M p ~ 7.2M Jup . As the stellar radius is R* = 6.7±0.3R Q JCarnev et al.ll2008bh . the true 
rotation velocity (on the equator) they derive is t> rot ~ 62 km s _1 . The orbital separation 
and eccentricity are a/ — 0.116 ± 0.01 AU and e = 0.25 ± 0.05, respectively. 

The initial angular momentum is practically that of the planet at its pre- RGB orbit, as 
the stellar angular momentum is negligible, 



J p0 = M p (GM* a 



|V2 



8.5 x 10 



M n 



M 



*0 



0.9M, 



1/2 



a 



2 AU") 



1/2 



g cm 2 s 1 , 



X 



,7.2M Jup , 

where M*q is the initial (pre-RGB) stellar mass, ao is the initial orbital separation, and we 
assume a circular pre-RGB orbit. 

The final angular momentum is carried by three components: The planet in its final 
(eccentric) orbit, J p ; the rotating envelope, J env ; and the mass that was expelled from the 
star, J w i n d. The present angular momentum of the planet is 



J p = M p [GMaj 



1/2 



1.8x10 



50 



M n 



7.2M Jup / V.0-8M 



'0 



1/2 



0.116 AU 



1/2 



g cm 2 s , 



wher e in the second equality we have substituted e = 0.25 and a = 0.116 AU (jSetiawan et al. 



12} 



2010). The current stellar mass is derived by taking a core mass of O.5M and an envelope 
m ass of Mhrctv ^ 0.3M ^. The envelope and co re masses are estimated based on the results 



of iDorman et al.l (I19931 ) and lD'cruz et al.l ( 119961 ). Although the envelope mass is estimated to 
be a little below O.3M , the large uncertainties concerning the angular momentum evolution 
justify using M H Benv = O.3M . 



The angular momentum of the rotating HB envelope is given by 

MHBenv 



J, 



HBenv 



R*v- 



rot 



1.7 x 10 



19 



0.3M m 



g cm 2 s , 



(3) 



where we took a = 0.01 from ISills fc Pinsonneaultl (120001) . a nd we have substituted R* 
6.7 R & and v rot = 62 km s" 1 as given by ISetiawan et al.l (120101 ) . 



To estimate the angular momentum carried by the wind we follow ISoker fc Harpaz 



(l2000l ) and assume that all the angular momentum was deposited on the RGB, and all 
mass-loss took place after the angular momentum was deposited. In the present scenario we 
propose, most of the mass-loss process took place while the planet was depositing its orbital 
angular momentum to the envelope. Such a process reduces the angular momentum carried 
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by the wind. On the other hand, the same scenario implies that as the planet deposited its 
angular momentum to the outer regions of the envelope, such that there was not enough 
time for convection to redistribute the angular momentum in the envelope. Such a process 
increased the angular mome ntum carried by t he win d. Over all there are large uncertainties, 



and we use the approach of iSoker fc Harpaa (120001 ). where more details are given. 



We assume the ratio of the HB envelope mass to that of the initial (before the interaction 
with the planet started) RGB envelope mass to be M H Benv/^RGBenv — 0.65. For example, 
this is the ratio for an RGB stellar mass of O.93M and a core mass of O.49M that give 
A^RGBenv = O.44M , and for a present envelope mass of MnBenv = O.28M . The total angular 
momentum of the envelope after deposition by the planet an d before mass-loss have started, 



based on the present angular momentum of the envelope, is (iSoker fc Harpaa l2000f ) 



-^HBenv \ ^ Q 1 n 50 ( MHBenvA ( ^HBenv M^RGBenv \ 6 2 -J 



j j ji£jciiv = 9 3 x 10 5U aacnv J1DU11V/ ^°™v „ cm ^ g-1 (A\ 

envR " em \M RGBcn J AdXiU \o.3MJ{ 0.65 J gCm ' W 

where the value of J cnv was taken from equation ([3]). The an gular momentum carried by the 



wind is J W md = ^cnvR — ^cnv The parameter 5 is derived by ISoker fc Harpaa (120001 ). which 
considered the range 4 < 5 < 7; we take here 5 = 6, but the uncertainties should be kept in 
mind. 

The total angular momentum carried by the different components after the planet 
started depositing its angular momentum to the envelope is not much below the estimated 
initial orbital angular momentum. 

Jcnv + ^wind + Jp ~ 4 X 10 5 ° g CU1 2 S^ 1 < J p0 . (5) 

Given the uncertainties, we can safely conclude that the total angular momentum of the 
system is about equal to the initial orbital angular momentum of the planet. This suggests 
that if a third body was present in the system its angular momentum was small relative to 
the initial angular momentum of the observed planet. Namely, it was a low mass planet 
and/or much closer to the star, and was swallowed earlier in the evolution. In any case, it 
did not play a dynamical role in the interaction between the star and the observed planet 
during the star transition from the RGB to the HB. This important conclusion will be used 
in section HI 



3. TIMESCALES CONSIDERATIONS 



We now show that a secular (regular) RGB evolution with a planet around it cannot 
lead to the present status of HIP 13044. The planet starts its journey towards the RGB 
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star when tidal interaction becomes strong enough to reduce the spiraling-in time below 
the remaining stellar evolution time on the RGB. For a planet mas s of M p ~ .01M < , n this 
occurs when the RGB radius reaches a value of -Rrgb — 0.25 — 0.4a (ISokerlll996l eq. 6 there; 
Villaver fc Livid l2009r iNordhaus et al.ll2010f ). After the spiral-in process has started, the 
process is accelerated tremendously. For a typical stellar parameters of a low mass star on 
the tip of the R GB, the spir alin g in time due to tidal in teraction when the planet is outside 
the envelope is (jSokerl Il996l and IVillaver fc Livid 12009k where the weak dependance on the 
other stellar parameters can be found) 



10 fc 



0.01M,c 



-i 



4-Rrgb 



yr. 



(6) 



When the planet reaches the RGB stellar surface the spiraling-in time is r s = r in (a = 
i? RGB ) ~ 10 yr. After the planet enters deep into the envelope, gravitational drag will 
accelerate the spiraling-in pr ocess, and the spiraling-in time becomes T; n (a < < 1 yr 
jNordhaus & Blackman1l2006h . 



For a planet to end at an orbital separation of aj = 0.116 AU the envelope should shrink 
to be less than cif within a timescale of t s ~ 10 yr. This is ~ 10 times the dynamical time 
scale of an RGB star, and shorter than the thermal time scale of the envelope r t h- e nv = 
GM RGB M env /(R RGB L ASGB ) ~ 100 yr, where R RGB = 0.5 AU and L A sgb = 1000L Q . 

It is unlikely that the interaction of the planet with the envelope by itself can cause 
the envelope to shrink over such a short time from ~ 0.5 AU to < 0.1 AU. If the core does 
not change the star stays a RGB star, that with an envelope mass of ~ O.3M has a radius 
much larger than 0.1 AU. The conclusion is that both the RGB core and envelope must 
be vigorously perturbed over a dynamical time scale while the planet is spiraling in. An 
interaction over a short time scale can account for the eccentric orbit of the planet in HIP 
13044 as well, as tidal and/or drag interaction over many orbits will circularize the orbit. 



4. DISCUSSION AND SUMMARY 

In section [2] we concluded, based on angular momentum considerations, that no tertiary 
object more massive than the planet was interacting with the RGB progenitor of HIP 13044, 
unless it had a much smaller orbital radius. In section [3] we argued that the interaction of 
the planet with the envelope must have been on a very short time scale of < 100 yr. We 
here present a speculative suggestion to account for these two conclusions. 

We speculate that the core helium flash caused the envelope of the RGB progenitor to 
expand by a factor of ~ 3 — 4 for a period of ~ 100 yr. In this scenario the planet was 
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orbiting the RGB progenitor at an orbital separation of a® ~ 2 — 4 AU. This speculative 
brief and large expansion of the RGB star during the core helium flash is composed of 
two steps. In the first step a small fraction of the energy released by the hydrogen that is 
ignited in the outer parts of the core, is transferred outside the core. This process is hard to 
study as it requires sophisticated 3D numerical study of the core helium flash. This energy 
deposition was n ot found until now in numerical simulations of core helium flash (e.g., 
Mocak et al.ll201ll ). Despite that, we here present some arguments that should motivate 
future studies to look for such an effect, in particular in rotating cores (that might have 
been spu n up by an inner pl anet early on the RGB). In that respect we note the comment 
made by iMocak et al.l (120081 ) that many processes following the core helium flash have some 
known inconsistency that indicate that the core helium flash is not fully understood. In the 
second step the energy deposited at the base of the envelope causes the envelope its large 
and brief expansion. Below we show that this is indeed the case. 

Along the RGB the star is powered by a hydrogen burning shell surrounding the almost 
pure helium core. When a temperature of little over 10 8 i^ is reached in the core, helium 
is ignited explosively. Because of neutrino cooling prior to the ignition, the ignition itself 
occurs off-center (e.g.. lMocakll2009l ). This core helium flash releases a vast amount of energy 
that cannot be carried radiatively, and thus convection is triggered in the core. It is thought 
that in most cases the He-burning and convective region reach up to the Hydrogen shell. 



becau se the H-burning shell provides an entropy barrier against mixing ( [Campbell et al. 



2010l and references therein). ICampbell et al.l (120101 ) point out that in solar- mass RGB stars 
of primordial or hyper-meta l-poor ( [Fe/H] < —5.0) , mixing might occur after all. There are 
two reasons for this mixing (ICampbell et al.ll2010l ). First, the core helium flash starts much 
farther away from the center in these low metallicity models than in solar metallicity stars. 
Second, the entropy barrier at the H-shell is much weaker in stars of very low metallicity 



becau se the H-burning shell almost switches off at this stage of evolution (IFujimoto et al. 



19901 ) . As a result of this mixing caused by the violent core helium flash in low metallicity 
stars (ICampbell fc Lattanziol 120081 ; ISuda fc Fujimotol I2OIO1 for \FelH] < - 2.5). ignition 



of large amount o f hydrogen occurs in these RGB stars (IMocak et al.l 120081 ; iMocakl 12009 
Mocak et allboioh . 



Mocak et al.l (120101 ) present a calculation of a core helium flash followed by hydrogen 
ignition. Over the first year the hydrogen burning provide ~lx 10 48 erg (see their fig. 1). 
After a year the hydrogen burning luminosit y is Lh ~ 1O 6 L . The hu ge energy prod uction 
by the hydrogen burning (IMocak et al.ll2010l ) and the core convection ([Blocker Ill999l ) decay 
over a time scale of ~ 10 — 100 years. Most of this energy stays in the core, and causes the 
core to swell. We now show that it is sufficient that ~ 5% — 10% of the energy released by 
the hydrogen burning leaks to the envelope to cause a substantial envelope expansion. 
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W e run a spherical evolutionary stellar code based on the one used by lHarpaz &: Kovetz 
( Il98ll ) with updated opacities. The initial composition used is X=0.689, Y=0.31, and 
Z=0.001. The code is spherical, and cannot follow the mixing of hydrogen to the core 
or the deposition of energy from the core to the envelope, as these processes are highly non- 
spherical. At the tip of the RGB, just as helium ignition starts, we manually add an anergy 
of 8.5 x 10 46 erg, at the bottom of the envelope, just above the hydroge n burning she l l. Thi s 
is ~ 7% of the energy released from the hydrogen burning reported by iMocak et al.l ( 120101 ) . 



The duration of the energy injection was 7 years at a power of L in = 1O 5 L . The initial 
model is presented in Figure HJ while the model at the end of the manually energy injection 
phase is shown in Figure [2j In Figure [3] we show the evolution of the outer radius and the 
outer boundary of the convective region. For the tidal interaction the outer boundary of the 
convective region is important. In this calculation the outer radius of the convective region 
increases by a factor of ~ 4. After ~ 100 yr the star shrinks back to its original radius. This 
calculation does not include mass loss, that can lead to further envelope contraction. 

The spiraling in of the planet with the properties of t he HIP 1304 4 system when the 
inflated envelope reaches the planet orbit is ~ 10 — 100 yr (jSokeri Il996l ) . The planet starts 
spiraling-in on a time scale about equal to that of the inflated envelope. 

We suggest that during the core helium flash of low metallicity RGB stars, a small 
fraction (few percents) of the energy liberated by the hydrogen burning, that takes place 
in the outer regions of the core, is transferred to the envelope. As we showed above, this 
causes the outer region of the envelope to substantially expand. A substantial increase in 
radiu s is found in some c alculations of shell helium fla shes (thermal pulses) in AGB stars 
(e.g. 



Schlattl et al.ll200ll ; iBoothrovd & Sackmannlll988f ). 



As the planet spirals-in in such an inflated envelope, it enha nces the ma ss loss rate by 
depositing gravitational energy and by spinning-up the envelope (jSokeri 12004 ). The calcula- 
tion of the mass loss is complicated and beyond the scope of this paper, but is expected to 
be of the order of ~ O.1M based on the properties of HIP 13044. The rest of the envelope, 
in our scenario, shrinks below the orbital separation of the spiraling-in planet before the 
planet manages to spiral-in below ~ 0.1 AU, and the spiraling-in ceases before the envelope 
is lost. This envelope contraction is caused by the changes in the core following the core he- 
lium flash, and it is expedited by the rapid mass loss caused by the spiraling-in planet. The 
rapidly changing envelope properties imply a rapidly varying tidal interaction, that instead 
of circularizing the envelope causes the eccentricity to increase. The overall evolution lasts 
for several dynamical time scales. This is possible because of the energy that is transferred 
from the core flash to the envelope over a very short time scale. This is crucial for our 
proposed scenario to work. 
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Initial model t=0 




20 40 60 80 
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Fig. 1. — The structure of the envelope just before the addition of energy at the base of the envelope 
(above the hydrogen burning shell). Shown are the mass (A/©), temperature (K) and density ( g cm -3 ), 
and the entropy (relative units). Convective regions are where the entropy decreases outward. 
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Fig. 2. — Like Figure 1, but after 7 years of energy injection at the base of the envelope at a power of 
L in = 1O 5 L . 
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Fig. 3. — The outer radius of the envelope (thin red line) and the outer boundary of the convective region 
(thick blue line) as function of time. Upper and lower panels show the same calculation but for different time 
spans. The noise in the graph (wiggling of the lines) demonstrates the numerical limitation of the code. It 
has a finite number of numerical shells, and the rapid increase in radius is done out of thermal equilibrium. 
These effects cause the noise. 
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We expect only low metallicity low-mass stars (Pop II) to experience the ~ 100 yr long 
inflated envelope phase. The total number of such objects in all globular clusters is expected 
to be < 1. Therefore, it will be extremely hard to find such stars in globular clusters and 
in the field. Even if found, they can be easily confused with AGB stars, unless they are 
followed for tens of years. In any case, the planet orbiting the red HB star HIP 13044 shows 
that planets can exist in globular clusters, and they can influence the evolution of the star. 
In particular they can increase the mass-loss rate and lead to the formation of blue HB stars 



(jSokerl Il998l ). We therefore suggest that globular clusters with a large population of blue 



HB stars be a prime target for exoplanet research. 

Another possible p rocess that might h ave occurred in this system is that it started as 



a multi-planet system (IBear &: Sokerll2010l ). There was an inner planet that was engulfed 
earlier on the RGB. This planet spiralled inward all the way to the core vicinity, and might 
have spun-up the core. Future 3D simulations of core helium flashes should check whether 
core rotation can facilitate the transfer of energy from the flashing core to the envelope. 

The Research was supported in part by the N. Haar and R. Zinn Research fund at the 
Technion, the Israel Science Foundation, and the Center for Absorption in Science, Ministry 
of Immigrant Absorption, State of Israel. 
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